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1. Introduction 

C/3 | Suppose f(x) is a real-analytic function denned in some neighborhood of the 

origin in R n . Consider the integral 
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a/(x V(^)^ (1-1) 



Here <f>{x) is a smooth bump function denned on a neighborhood of the origin and A is 
a real parameter whose absolute value we assume to be large. If V/(0) 7^ 0, then by 
repeated integrations by parts (see [S] Ch 8 for example), for any iV one has an estimate 
1 1\ I < C/^jv I A| for appropriate constants C/^jv- In the case where V/(0) = 0, that 
is, when / has a critical point at the origin, it can be proven (see [M]) using Hironaka's 
resolution of singularities [H1]-[H2] that if the support of <fi is contained in a sufficiently 
small neighborhood of the origin, then as A —>■ 00, I\ has an asymptotic expansion of the 
^ ' form 
(N 



A/(0) 

a i=0 

Here the sum in a goes over an increasing arithmetic progression of positive rational 
numbers, and the aj )Q are distributions with respect to the cutoff <j). We refer to the 
excellent resource [AGV] for further results along these lines. In Theorem 1.2, we will 
provide another proof of the existence of the expansion (1.2) using an elementary resolution 
of singularities theorem deriving from [Gr], thereby avoiding the use of Hironaka's theorem 
or other nonelementary techniques. Illustrating this elementary method of proving (1.2) 
and related results while giving the precise estimates of Theorems 1.1-1.3 can be viewed as 
the main purpose of this paper. In Theorem 1.4, we will also give an elementary proof of the 
well-known result of Atiyah [At] and Bernstein-Gelfand [BGe] concerning the meromorphy 
of integrals of f z for nonnegative real-analytic /. It should be pointed out that there exist 
other elementary resolution of singularities algorithms which have been used for various 
purposes, notably [BiMi]. 

The analysis of (1.1) is closely related to the analysis of sublevel set integrals (see 
Ch 6-7 of [AGV] for example). Namely, let f(x) be as above and assume /(0) = 0. We 
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consider the integrals 

J t = (f>(x) dx (1.3) 

J{x:0<f(x)<t} 

Theorem 1.1 below, also proved using an extension of [Gr], will show that for small t > 
we have an asymptotic expansion 

n-l 

a i=0 

Furthermore, each ^q^t will be seen to have asymptotic expansion given by termwise 
differentiation of (1.4). This appears to only have been explicitly done in the case where 
/ has an isolated zero at the origin [J] [Va] , although one should note that if one is willing 
to assume Hironaka's results it can be proved without a huge amount of difficulty. [L] 
considers some related issues in this subject. Once one knows the expansion (1.4), one 
can obtain the expansion (1.2) using well-known techniques. In the integral (1.1), after 
factoring out a e tX ^ ^ if necessary, one may work under the assumption that /(0) = 0. 
In this situation one first integrates over a given level set f = t and then integrates with 
respect to t. Consequently, the integral may be rewritten as 

Here J t denotes the analogue of J t with / replaced by — /, and 7(£) denotes a bump 
function such that r )(f{x)) = 1 for all x G supp{4>). With a little care, one can then 
substitute (1.4) for / and — / respectively into (1.5), integrate term by term, and obtain 
(1.2). We will do this rigorously in the proof of Theorem 1.2 below. 

It turns out that there are some natural generalizations of (1.4) which are not 
any harder to prove using resolution of singularities than (1.4) itself. For example suppose 
A = {x E R n : gi(x) > 0, gk(x) > 0} where the g t are real-analytic. Assume further 
that G bd(A). Then one can obtain an asymptotic expansion for /| xeJ 4.o</( x )<t} ( P( X ) dx 
in the same fashion that one obtains (1.4). Furthermore, if one has several real-analytic 
functions fi(x), fi(x) each satisfying fi(0) = 0, then one can similarly obtain an asymp- 
totic expansion for f^ xeA . 0< f 1 ^ o<fi(x)<t} ^( x ) dx - However, the issue of exhibiting an 
asymptotic expansion for J {xeA:0<fl{x)<tl: ^^ 0<fl{x)<tl} 4>(x) dx in appears to be 

quite a bit harder, and will not be addressed here. Our first result is the following, to be 
proven in section 3 using the elementary resolution of singularities theorem of section 2. 

Theorem 1.1: Suppose fi(x), fi(x) and gi(x), gk{x) are real-analytic functions 
defined on a neighborhood of the origin in R n , with fi(0) = for all i. Let A = 
{x G R n : gi(x) > 0, ...,gk(x) > 0}, and assume G bd(A). There is a neighbor- 
hood V of the origin such that if <f>(x) is a C°° function supported in V, then J t = 
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I{xeA:0<M x )<t, ... ,o<Mx)<t} <K X ) dx has an asymptotic expansion given by 

n-l 

Jt = E ^,«(0)* Q ln(t) J + £ a (t) (1-6) 

Here the a range over an arithmetic progression of positive rational numbers depending 
on the fi and the g^. Let Z = {x G cl(A) : fi(x) = for some i.} There are M > and 
A a > (depending on the fi and Oj) such that each bi jCe is a distribution with respect to 
4>, supported on Z and satisfying 

\bi, a ((j))\ < A a sup sup (1.7a) 

\P\<Ma xeZ 

The error term E a (t) is such that there is e > and C a > such that if < m < a, then 

lm 

\ E a (t)\<C a sup sup 1^0(^)1^+— (1.76) 

at \P\<Ma x 



The theorem of this paper regarding oscillatory integrals is as follows: 

Theorem 1.2: Suppose f(x) is a real-analytic function defined on a neighborhood of the 
origin in R n with /(0) = 0, and let A and F be as in Theorem 1.1. Then as A — > oo, 

/a = £4 e lX f( x }(f)(x) dx has an asymptotic expansion 

n-l 

^^a l)a (0)A- Q ln(A) l + K(A) (1.8) 

a<a i=0 

The a range over an arithmetic progression of positive rational numbers depending on / 
and the gj. Let Z = {x G cl(A) : f(x) = 0}. There are M' > and > such that each 
cii ta is a distribution with respect to 0, supported on Z and satisfying 

k Q (0)|<4, sup sup|^0(x)| (1.9a) 

\P\<M'a xEZ 

The error term E' a (t) is such that there is e' > and C' m a > such that for any m one 
has 

tKWKC sup sup \d^{x)\\-*-'-™ (1.96) 

« A |/3| <M'(a+m) a; 

Proof: As in (1.5) we write 
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^e at 7 (t)dt+ / ?p. e - iXt >y(t)dt (1.10) 
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As before J t denotes the analogue of J t with / replaced by — / and 7(f) has compact 
support and is equal to 1 on a neighborhood of the origin. The two terms of (1.10) are 
done in the same fashion, so we focus our attention on the first term. By Theorem 1.1, for 
a given k we can write 

a <k-l i=0 

Here the Bi^ a are obtained from the bi yCe by performing the appropriate term-by-term 
differentiation. In order to be able to differentiate the expansion (1.8) up to A~ a a total 
of m times, m > 0, we need to insert the above expansion for k = a + m + 2 into the term 
of (1.10). We get 

n — l „oo poo irp 

h= £ J2 B ^ t a Ht) i e iXt 1 {t)dt+ "+ m+2 (t)e iXt >y(t) dt (1.11) 

Q < a+m+ll=0 Jo Jo at 

It is well-known (see [F]) that for any / > one has 

/ e lXt t a \n(t) m 1 (t)dt= . \/T > +0(\- 1 ) (1.12) 

As a result, for any /, (1.11) becomes 

'*= E E^-wli^+f to( t )^ 7W * +0 (A-) (i.i3) 

a<a+m+l i=0 V ' " y> 

Equation (1.13) will give the desired expression for the mth derivative of I\S expansion 
up to order A _a . For in (1.12) the 0(X~ l ) behaves as needed under differentiation (again 
see [F]). As for the d ^ a + m+2 term in (1.13), one can differentiate J °° dEa %™ +2 (t)e iXt -f(t) dt 
under the integral sign, obtaining a power of it for each of the m A-derivatives taken. Then 
one does a + m + 1 integrations by parts in t, integrating the e lXt factor and differentiating 
the rest. Equation (1.76) ensures the left endpoint terms disappear. We obtain 

, d m f°° dEg+rn+2 iXt -g-m-1 _ | d a+m+2 E a+m+2 , , , 

l dX^J Q ~^^ {t)e ^)^|<C a , m A sup I dfa+m+2 (t)\ 

<C' a ^ sup sup \d^{x)\\- a - m - x (1.14) 

\/3\<M'(a+m) x 

The last inequality follows from the ,, 0+m a f t a m+2 (t) case of (1.76). Inserting this back 

into (1.13) gives the desired result. This completes the proof of Theorem 1.2. 

Next, we focus on the situation where / > on A. Where r is now a large 
positive parameter, we consider the Laplace Transform-like object defined by 

e- Tf(x) (P(x)dx (1.15) 
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We have the following theorem regarding L T . 



Theorem 1.3: Suppose f(x) is a real-analytic function defined on a neighborhood of the 
origin in R n with /(0) = 0, and let A, V, and Z be as in Theorem 1.1. Suppose on a 
sufficiently small neighborhood U of the origin we have that f(x) > on U D A. Then if 
the support of 4> is contained in V fl U, as r — > oo, L T has an asymptotic expansion 



n-l 



L - = E E <M^~ a Wr) 1 + (1.16) 



a<a i=0 



The a range over an arithmetic progression of positive rational numbers depending on / 
and the gj. There are M" > and A n a > such that each di jCe is a distribution with 
respect to 0, supported on Z and satisfying 



Ka(0)l<< SUp SUp|^0(x)| (1.17) 

|/3|<M"a x£Z 

The error term E"(t) is such that there is e" > and > such that for any m one 
has 

Jm 

tKWK^, a sup sup |^0(x)| r— (1.18) 

°" |/3|<M"(a+m) x 



Proof: The proof is essentially a repeat of that of Theorem 1.2. Namely, we have 

/•OO o T 

Lr = i, l£ e ~ Ttl{t)dt (L19) 

Because of the / > condition, there is only one term instead of two this time. Instead of 
using (1.12), here we use 

f°° 8 i Yin 4- \ \ 

I e-^Hty 1 {t)dt = — i ^ T ±+0{r l ) (1.20) 

Otherwise the proof is identical to that of Theorem 1.2 so we omit the details. Incidentally, 
equation (1.20) is somewhat easier to prove than (1.12). When i = one can write (1.20) 

as as 

/•OO /-OO 

/ e~ Tt t a dt- / e- Tt t a (l--f(t))dt (1.21) 
Jo Jo 

The first term is exactly r ^_ZtP ■ On the other hand, via repeated integrations by parts, 
the second term and its r derivatives are seen to be 0{r l ) for any /. The % > case follows 
from differentiating (1.21) under the integral with respect to a. 

Theorem 1.1 also gives as a relatively straightforward consequence the following 
theorem of Atiyah [At]; a similar result is due to Bernstein-Gelfand [BGe]: 
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Theorem 1.4: Let /, A, U, and V be as in Theorem 1.3. Define F^{z) = J A f(x) z (j)(x) dx. 
Then if the support of <fi is contained in V D U, the function F^z), initially defined as a 
holomorphic function of z on Re(z) > 0, extends to a meromorphic function on all of C. 
The poles of this extension are located on an arithmetic progression of negative rational 
numbers depending on / and the gj, and each pole is of order at most the dimension n. 

Proof: Replacing / by cf for an appropriate c > if necessary, we may assume that 
\f \ < 1 on supp(4>). Analogous to (1.19), we have 

r 1 d T 

F 4> (z)= J t z ^dt (1.22) 
Inserting (1.11) into (1.22), one obtains 

^) = EEW) [ 1 r +z Htydt+ f^tydt (i.23) 

a<ai=0 J ° J ° M 

When Re(z) > —1, each of the integrations of the sum in (1.23) can be performed directly, 
and we obtain 

n dV 

a<a i=l J ° 

Note that the sum in (1.24) automatically extends to a meromorphic function on C with 
poles of order at most n. As for the error term, one can rewrite it as f Q F a (t)t z+a ~ 1+e dt, 
where in view of (1.76) F a (t) is bounded. Hence on Re(z) > —a + 1 — e, the error term is 
an average of locally uniformly bounded analytic functions. Hence it is itself an analytic 
function on Re(z) > —a + l — e. Since a can be made arbitrarily large, the theorem follows. 

The remainder of the paper is organized as follows. In section 2, a version of 
the elementary resolution of singularities algorithm of [Gr] will be developed that will be 
appropriate for proving the type of theorems of this paper. It may also be useful for other 
purposes as well. Section 3 will be devoted to proving Theorem 1.1 using the algorithm of 
section 2. 

2. A resolution of singularities theorem 

We now give some terminology from [Gr]: 

Definition: We say that a function g : A C R n — > R n a quasitranslation if there is a real 
analytic function r(x) of n — 1 variables such that g(x) = (g\(x), g n (x)), where for some 
j we have gj(x) = Xj — r(xi, ...Xj-i, Xj + \, ...,x n ) and where gi(x) = Xi for all % ^ j. In 
other words g(x) is a translation in the Xj variable when the others are fixed. 
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Definition: We call a function m : A C R n — > R n an invertible monomial map if there 
are nonnegative integers {«ij}" J=1 such that the matrix (cty) is invertible and m(x) = 
(mi(x), ...,m n (x)) where rrii(x) = x^ 1 ....x% in . The matrix (ct^) being invertible ensures 
that ft is a bijection on {x : x\ > for all /}. 

Definition: Let E = {x : Xi > for all i}. If h(x) is a bounded, nonnegative, compactly 
supported function on E, we say h(x) is a quasibump function if is of the following 
form: 

/»(*) = a(s)nMc«(s)^) (2.1) 

Here pi(x), qi(x) are monomials, a(x) G C°°(cl(E)), the are nonvanishing real-analytic 
functions defined on a neighborhood of supp(h), and are nonnegative functions in 

C°°(R) such that there are constants c\ > Co > with each 6/ (x) = 1 for x < Co and 
6/(x) = for x > c\. 

The main theorem from [Gr] is as follows: 

Main Theorem of [Gr]: Let f(x) be a real-analytic function defined in a neighborhood of 
the origin in R n . Then there is a neighborhood U of the origin such that if <f>(x) G C£°(U) 
is nonnegative with 0(0) > 0, then can be written (up to a set of measure zero) as a 
finite sum £\ </>i(x) of nonnegative functions such that for all z, G supp((pi) and supp(<pi) 
is a subset of one of the 2 n closed quadrants defined by the hyperplanes {x m = 0}. The 
following properties hold. 

(1) For each % there are bounded open sets D^,...,D^ i , and maps g},..., g^, each a reflection, 
translation, invertible monomial map, or quasitranslation, such that D® = {x : (pi(x) > 0} 
and such that each gj is a real-analytic diffeomorphism from Df to D\~ . The function gj 
extends to a neie; hborhood N- of cZ(Z^) with g\ (N?) C N- for j > 1 and ^(^Z) C t7. 

(2) Let E = {x : x t > for all i} and = ^ o .... o gj?\ Then C -E, and there is a 
quasibump function $j such that x D k % { x ){4>i ° ^i( x )) = &i( x )- 

(3) G AT^ with *i(0) = 0. 

(4) On N** , the functions / o det(^i), and each jth component function is of the 
form c(x)m(x), where m(x) is a monomial and c(x) is nonvanishing. 

The corollary to the main theorem of [Gr] says that one can resolve several func- 
tions simultaneously in such a way that the resolution satisfies the conclusions of the main 
theorem. However, these theorems are not precisely what is needed for the arguments of 
this paper, because here a quasibump function is not the appropriate form for the function 
of part (2). Instead we will need the following: 

Theorem 2.1: If in the main theorem of [Gr] and its corollary, if one replaces the condition 
in part (2) that $>i(x) = x D k t {.%){4>i°^ i( x )) is a quasibump function with the condition that 
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for some rectangle Ri = (0, a\) x ... x (0, a l n ) the function &i(x) is of the form XRi ( x )l( x )i 
where 7 (a?) is a C°° function on cl(Ri), then the rest of the main theorem and its corollary 
respectively still holds. 

One way of looking at Theorem 2.1 is that in the blown- up coordinates, one can replace 
the issues coming from the singularities of the quasibump function by the issues coming 
from the jumps in the characteristic function of Ri. These latter issues will turn out to 
cause no problems in the analysis of the integral quantities of this paper. Theorem 2.1 will 
be a consequence of the following, which we will prove later in this section. 

Theorem 2.2: Let <&(x) be a quasibump function, and let pi(x) and qi(x) be as in the 
definition (2.1) of quasibump function applied to $>(x). Then there is a J such that for 
1 < j < J and 1 < k < n there are invertible monomial maps 9j(x), quasibump functions 
Qj(x) of the form ]T^i a (f^) and sets Fj = {x e E : < 1 for 1 < k < n} where 

rjk(x) and Sjk(x) are monomials, such that the following hold. 

(1) Up to a set of measure zero one has a decomposition 

j 

<S>(x) = J2®(x)Qj(x)XF j (x) (2-2) 
i=i 

(2) Each ^4#TT and each Pl \ 9 ^ x ^. is a monomial. 

(3) For each j there is some rectangle Rj = (0, cr}) x ... x (0, a J n ) such that XFj(gj{x)) = 

Comment: Note that <&(g(x)) = a(gj(x)) YliLi h (q ( 9j (x) ) ^ ^ J j ) , a and c\ smooth, and 

that Qj(gj(x)) = nl^i a ( q j -l(l (*c)) i )• Hence (2) implies that §(gj(x)) and each Qj(gj(x)) 
are smooth. Thus by (3) , for each j there is and a smooth 7j (x) on Rj with 

^(gj(x))Qj(gj(x))x Fj (gj(x)) = XR j {x)lj{x) (2.3) 
Proof that Theorem 2.2 implies Theorem 2.1: 

Suppose Theorem 2.2 is known to hold, and let f(x) and U be as in the main theorem 
(or its corollary). Suppose 4>{x) is a bump function defined in U, and let 4> = J2i4>i be 
the decomposition of the main theorem (or its corollary). Let <&i(x) be the quasibump 
function and the composition of coordinate changes associated to 4>i- Let Qij( x ), Fij, 
and gij be as given by Theorem 2.2 applied to $>i(x). I claim that the decomposition 
4>(x) = J2ij 4>i( x )Qij('$~ 1 (x))xF ij ('$i 1 (x)) satisfies the conditions of the main theorem 
of [Gr], with coordinate changes gj, g i i , g^, if one specifies the domains Di and Ni as 
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follows. At the kith level, the domain is defined to be n F^ n {x : Qij(x) > 0}. 
For k < ki, the domains are sucessively defined by F>\~ x = g^{D^-). At the final ki + 1th 
level, one puts D^j +1 = g^D^. For j < ki one can define to just be N? , and then 
to be ^JVj. 

With the above definitions, part (1) of the main theorem is readily seen to hold. 
As for (2), since x D k i ( x )(4>i ° ^i( x )) = ®i(x), we have 

i 

XflJiW^i o ^i{x))Qi j {x) X F ij {x) = ^i{x)Qi j {x) X F ij (x) (2.4) 

Therefore, 

^■n{x:Q l3 (x)>o}nF i3 ( 1 )^' ° ViWQijWXFuix) = ^i{x)Qi j (x) X F ij {x) 
Equvialently, 

X D ki{x){(j)i o ^ i (x))Q ij (x)xF ij (x) = ^i(x)Q ij (x)xF ij (x) 

Composing with gij, we have 

X D k i+ i(x)((j)i o o g lJ (x))Q tJ (g tJ (x))xF iJ (gtj(x)) = <S> l (g l j(x))Q lJ (g lJ (x))xF lj (gtj(x)) 

(2.5) 

If one lets = o ^ and (pij(x) = (pi(x)Qij($ \ 1 (x))xF i:) (^i 1 (x)), from (2.5) and the 
assumption (2.3), for an appropriate rectangle Rij one gets 

X^+iix^iji^ijix)) = x D k i+ i(x)((j)i o ^i o g ij (x))Qi j (gi j (x))xF ij (gij(x)) 

ij ij 

= (g^ (x))Qij (g tJ (x))x Flj {gij (x) ) 

= XR ij (xHj(x) (2.6) 

This gives the version of part (2) of the main theorem that is needed in Theorem 2.1. Parts 
(3) and (4) are immediate, and we are done. 

Proof of Theorem 2.2: The proof is by induction on the dimension n. When n = 1, 
since 1-dimensional quasibump functions are smooth already the proof is straightfoward. 
Namely, one selects T\ such that supp^) C [0, r{\. We then let there be a single Fj = 
X(o,n)( x ) an d a single Qj(x) = a(x) where a(x) = 1 on [0,ri]. The corresponding gj is 
just the identity map, and the case n = 1 follows. Assume now we that know Theorem 
2.2 in dimension n — 1 and the hypotheses of Theorem 2.2 hold for some n-dimensional 
situation. We break into two cases. 



Case 1: The first case is when there are distinct monomials t(x) and u(x) and positive 
constants c\ and C2 such that whenever <&(x) ^ we have 

t(x) 

ci < 44 < c 2 2.7 
u(x) 
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Write t(x) = Yliei^i an< ^ u ( x ) = Yliei' x i *' where each /; > 0. We can assume that 
I n I' = and I U I' 7^ 0. We change variables as follows. For a given z e / U let 
m i = rij6/u/',j/i ^' anc * let x * = ^ mi - If * ^ -f U let = Let x = <?i(y) be this 
coordinate change, and let / denote Ilie/Li/' ^- Then we have 

*G/i(y)) = II^ «G/i(y)) = II^ ( 2 - 8 ) 

iei iei' 

Consequently, if t 1 (y) = Yl ieI yi = [t(gi(y))]i , and = Fliej' 2/* = Hdiiv))}^ , equa- 

tion (2.7) says that whenever &(gi(y)) ^ one has 

c t < My) < i (2 . 9) 

wi(y) 

Effectively, we have reduced to the case when each /; is 1. We now prove Theorem 2.2 for 
case 1 by induction on m = min(|I|, We start with the case where m = 0. In this case, 
either u\{y) or t\(y) is a nonconstant monomial and by (2.9) whenever <&(gi(y)) ^ that 
monomial is bounded below. Since the support of a quasibump function is also bounded, 
every yi appearing in this monomial is therefore bounded below on the support 3>(gi(y)). 
Specifically, there is some yi and some constant c such that yi > c whenever &(gi(y)) is 
nonzero. 

Since <&(x) is a quasibump function, so is 3>(gi(y)). As in (2.1), we write 3>(gi(y)) as 



k 



*{y)J[Uci{y) p ^) (2.10) 

Because yi is bounded below, one can incorporate any powers of yi appearing in each pi(y) 
and qi(y) into the associated ci(y). Thus we may assume that the pi(y) and qi(y) do not 
depend on j/j. Let c > be a constant such that nf=i ^( c f7[fy) = 1 on the support of 
<3>(<7i(j/)). So we have 



*(gi(y)) = HgM)\[Hc p -^) ( 2 - n ) 

We apply the induction hypothesis to nf=i M c f7(fy)' a function of n — 1 variables. Let 
Qj and Fj be as in Theorem 2.2 applied to this function, and let gj(z) be the associated 
invertible monomial map, called gj(x) in the statement of that theorem. Then 

nw^)=E(nMc^)Q,(^w) (2.12) 

Combining the last two equations gives 

*a/i(y)) = E(^^(v))n 6 '( c ?M)]^^(i/)) ( 2 - 13 ) 
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Using this again on the bracketed expression, we have 

®(gi(y)) = £*toi(v))Qj(v)xF,(y) (2.14) 

3 

This will lead to our needed decomposition for ^(g 1 (y)). Next, note that we have 

*(giog*(z)) =Y J H9i{9%z)))Q j {g%z)) XFj {g 2 j {z)) (2.15) 

3 

Since Theorem 2.2 holds for YlLi b i( c ^m) and 9j(z), each El Wr^ in $ (fi(f , (^))) is a 
monomial. Similarly, if as in Theorem 2.2 we write Qj(y) = Yli bji(cji(y) g J )(yj )> we a l so 
have that each Pj ^ g j^^ appearing in Qj(g J 2 (z)) is a monomial. So since Xf,- (<72(<z)) = 
XRj(z), equation (2.15) shows that Theorem 2.2 holds for $, with the associated invert- 
ible monomial maps given by g\ o g 3 2 . Hence we are done with the case where m = 
min(\I\, \I'\) = 0. 

Now we assume min(\I\, \I'\) = m > 0, and that we know the theorem for 
min(\I\, = m — 1. For this fixed m, we induct on max(\I\, \I'\) which we denote by 
M. The initial step M = m is done the same way as the inductive step, so we assume 
either M = m or that M > m and that we know the result for M — 1. Let %\ e / and 
let %2 € i 7 . Let ct(t) be a C°° function on [0, oo) that is equal to 1 for small enough t and 
which satisfies a(t) + a{\) = 1. In particular, 

«(— ) + «(— ) = 1 ( 2 - 16 ) 

yi 2 v%i 

Correspondingly, we decompose <&{g\(y)) = $>i(y) + <&2(j/), where 

$i(y) = <f>(gi(y))*(^), *M = Hgi(y)M-) (2.17) 

Let G\(y) be the invertible monomial map whose iith component is yi x yi 2 , and whose ith 
component is yi for all % 7^ %\. Similarly, let Gi{\j) be the invertible monomial map whose 
Z2th component is y^y^, and whose ith component is yi for all % 7^ z 2 . We have 

SiCGifc/)) = ^feod^Ja^), * 2 (G 2 (y)) = <&(<7i o G 2 (y))a(y ia ) (2.18) 

Both $i(Gi(y)) and $2(^2(2/)) are quasibump functions for the following reasons. First, 
since <&(gi(y)) has compact support and for k = 1,2 the coordinate change Gk is in the 
yi k variable only, the &k(Gk(y)) have compact support in the yi variable for i ^ ik- When 
i = ik, the a(yi k ) factor ensures that 3>k(Gk(y)) has compact support in the yi h variable 
as well. We conclude that both &k(Gk(y)) are compactly supported. Furthermore, since 
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the Gk are invertible monomial maps and $ is a quasibump function, the <&k(Gk(y)) are 
automatically of the form (2.1). We conclude they are quasibump functions. 



Next, we translate the condition (2.9) in the new variables. Since appears in 
ti(y) and yi 2 appears in ui(y), both of degree 1, the expression for ^(^/(y)) is obtained 
from the expression for by removing yi 2 from the denominator, while ^(G^y)) ls 

obtained from the expression for by removing from the numerator, In either 

case, either m = min(\I\, or M = max( |/|, \I'\) is reduced by 1. So by the inductive 
hypothesis Theorem 2.2 applies to both &i(Gi(y)) and & 2 (G 2 (y)). Let the decompositions 
from Theorem 2.2 applied to these functions be given by 

$ 1 (G 1 (y)) = Y^®i(G 1 (y))Ql j (y) X Fl.(y), **(G 2 (y)) = ^* 2 (G 2 (y))Cg j (y)x F *(y) 
3 3 

(2.19) 

Pulling back by G\ and G2 respectively, we have 

^i(y) = J2^QU G i 1 y)XFliG^y\ ^ 2 (y) = J2^2(y)QUG2 1 y)XF^ 1 y) 

3 3 

(2.20) 

Since &(gi(y)) = &i(y) + $2(2/), we may add these, obtaining 

$(si(y)) = E $ i(^( G rWi^^^ (2.21) 

3 3 

Going back to the original x coordinates 

^x) = J2^i(9i\x))Q} J ((g l oG 1 )- 1 x)x F ii(9ioG 1 )- 1 x) 
3 

+ J2^2(9i 1 (x))Ql((g 1 oG 2 )- 1 x)xFi j ((9ioG 2 r 1 x) (2.22) 

3 

Equation (2.22) gives the required decomposition of the form (2.2) for $(x). If gjp gfj de- 
note the invertible monomial maps coming from Theorem 2.2 in (2.19), then the invertible 
monomial maps corresponding to the decomposition (2.22) are given by g\ o G\ o g\- or 

9l oG 2 o g%. Since ^\x)) = a (gSg^)*(s) and <5> 2 (g^(x)) = a (gSg^)*(x), 

the quasibump functions for (2.22) (corresponding to the Qj(x) in (2.2)) are given by 

a( ( ^ x \^ )Q}M gi o G 1 )~ 1 x) and a( { ( 9 U ( , X ^ 2 )QMgi ° G^x), while the sets corre- 

sponding to Fj in (2.2) are given by (#1 oGi)~ 1 Fl j or (gi oG 2 )~ 1 F? j . That (2.22) satisfies 
the conclusions of Theorem 2.2 is a direct consequence of the fact that the decomposition 
(2.19) does. This completes the proof in Case 1. 
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Case 2: The second case is when no expression of the form (2.7) holds. Since is a 
quasibump function, as in (2.1) we can write 



^x)=a(x)f[b l (c l (x)^-) (2.23) 
i=i qiy > 

Since $>(x) is compactly supported, we can multiply (2.23) through by []^ =1 a(x m ) for 
an appropriate function a(x) and not change the result. Hence we can assume the x m are 
amongst the j^jy- Write = nlLi 3 ^*' - Here an can be positive, negative, or zero. 
By definition of quasibump function, there is a constant c\ such that <&(x) is supported 
on {i G £ : < ci for all /}. If we define (j/i, y n ) coordinates by yi = ln(xi), and 

write the associated coordinate change as x = e{y), then $(e(y)) is supported on the set 
A given by 

n 

^ = r? l=1 {y G R n : J2 a ^ < ln (ci)} (2.24) 

i=i 

Define the set A' by 

n 

^ = nf =1 {2/ G R n : < °) ( 2 - 25 ) 

i=l 

Lemma: A' is nonempty. 

Proof: We will show that if A' = 0, then one must be in case 1 of this proof. Suppose 
A' = 0. Note that we may assume ln(ci) > 0; otherwise A' D A ^ 0. Since A' = 0, we 
have A = A — A' and therefore 

n 

A C U J l=1 {y G R n : < ^ < ln(ci)} (2.26) 

i=l 

I claim that there is some M > such that 

n 

A C {y G R n : -M < J2 a *iV* < Hci)} (2.27) 

i=i 

For suppose not. Then any M > we can find a, y l E A such that $^l=i c^z^ < — M. If y 
denotes the average - Xw=i y l , then by the convexity of A, y G A. Furthermore, since for 
all l\ and Z2 we have Y17=i a ihVi < l n ( c i)? taking this average leads to 

El \ - \ - , j — 1 M 

i=l ^ 2=1 i=l ^ ^ 

Hence if M is large enough, one has 

n 

Y j ailVi<® (2-28) 
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Hence y G A', contradicting that A' = 0. Therefore (2.27) must hold. This automatically 
implies that we are in case 1; for the y coordinates, (2.7) translates to the existence of 
numbers e±, e2, di,...,d n such that e\ < Y^i=i^iVi < e2 whenever (yi,...,y n ) G A. Thus 
(2.27) implies we in case 1 and the lemma is proved. 

Lemma: There is a vector v such that A C A' + v 

Proof: let w be any vector in A' . Then we have Yl7=i a u w i < f° r every I- If V £ -4) for 
any t > 0, the vector y + tw satisfies 

n n n n 

^2 a n(y + tw )i = ^2 a nVi + t(%2 auwi) < ln(ci) + t(y] auwi) 

Thus if t is large enough, one has Yl7=i a n(y + t w )i < for all y G A. In other words, 
y + tw G A'. Hence A C A' — tw and we are done. 

Since the coordinate functions x m are amongst the (see below (2.23)), the 
n inequalities {y m < 0} are amongst the defining inequalities for A' and therefore 

A' C {y : y m < for all m} (2.29) 

We now "triangulate" A'. To be precise, by (2.29) we see that A' C\ {y : Ylm=i y™ = ~ 1} 
is a bounded convex polyhedron and therefore up to a set of measure zero can be written 
as a finite union DjSj of n — 1-dimensional simplices. For a given j, we let Tj be the union 
of all lines containing the origin and passing through Sj. We then have 

A' = UjTj 

Since A C A' + v for an appropriate vector v, we have 

A C LI, (Tj + v) (2.30) 
Furthermore, each Tj can be written in the form 

n 

= nr =1 {y G R n : ^^yl/i < °} ( 2 - 31 ) 

i=i 

Equivalently, 

n n 

Tj + v = Df =1 {y G R n : h HjVi < Yl hil ^ 

i=i i=i 

We write rjij = Y^=i huj v ii so that the above becomes 

n 

Tj+v = nf =1 {y G R n : ^ fc^y* < 7/y} (2.32) 

i=l 
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Next, we move everything back to the original x coordinates. (Recall yi = ln(xi) for each 
i). Denote the associated coordinate change by x = e(y). Writing E = {x : Xi > for all 
z}, we have 

n 

e(Tj + v) = n? =l {x G E : ]Jx^ lj < exp(rjij)} (2.33) 

i=l 

By (2.30), e(A) C Uje(Tj +v), while in (2.24) we defined A so that the original quasibump 
function $(x) is nonzero only on e(A). Hence if we denote e(Tj + v) by Fj, we have 

For each j let Qj(x) be some bump function of the form Y\7=i a ( x i) that is equal to 1 on 
the support of $>(x). One then has 

$(x) = "£Hx)Qi(x)XF j (x) (2-34) 

j 

This will be the decomposition needed for Theorem 2.2. For a given j, the associated 
coordinate changes gj(z) are defined as follows. Let Hj be the matrix whose il entry 
is hnj, and let Mj = {edij}di=i a m atrix such that integer coordinates such that 
MjHj is iV times the identity matrix for a large integer N. If one does the substitution 

Xi = Yld=i z d dl3 1 then nr=i x f llJ becomes zj^ . We define gj(z) to be the map 

n n 
d=l d=l 

Note that by (2.33) we have 

g-'F, = (0, expCf)) x .... x (0, exp(^)) (2.35) 

Claim: The Fj and gj(z) defined this way satisfy the conclusions of Theorem 2.2. 

Proof: First, we check the gj(z) are invertible monomial maps; we have not shown that 
each €ijk is nonnegative. Since Fj C supp(<&), Fj has compact support. In particular each 
Xi is bounded on Fj. Pulling back to the z coordinates, this means each Y[ d =i z d d ^ is 
bounded on (0, exp(^jf)) x .... x (0, exp(^f)). This fact forces each eaij to be nonnegative. 

For on the curve z = {t ai , t Q "), ad > 0, the ith component of gj{z) is equal to ad£di J _ 
If some edij were negative, by choosing ad to be far larger than the other the image of 
the curve under gj would go off to infinity. We conclude that each €dij is nonnegative and 
therefore that gj(z) is an invertible monomial map. 

Next, observe that by (2.34), part (1) of Theorem 2.2 holds. As for (2), we must 
show that if pt and qt are as in (2.23), then each Vk \ 93 \ z Jl is a monomial. Recall that by 
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(2.31), we have Tj = n? =1 {y G R n : J27=i h UjVi < °} 5 and that each equation f^Jfyj < 1 
is one of the defining equations for A. Translating the statement that Tj C A into the x 
coordinates gives: For fixed j, if Y\7=i x i" j < 1 for a11 1 then < 1 for a11 k - Pulling 

back by g 2 , we have that if z G (0, exp{^-)) x .... x (0, exp(^)), then gjg^jj < 1 for all 
fc. Note that ^[g^]) is of the form rid=i f° r integers 5dkj- So exactly as in the last 
paragraph, we must have that 5dkj > for each d and k. Hence for each k, ^7^4^ is a 
monomial. Since j was arbitrary, we have part (2) of Theorem 2.2. Equation (2.35) gives 
part (3) and we are done with the proof of Theorem 2.2. 

3. Proof of Theorem 1.1 

Most of this section will be devoted to proving Theorem 1.1. We will use the 
resolution of singularities algorithm as described in Theorem 2.1 to reduce consideration 
to the simpler situation where there is only one fi(x), and where that fi(x) is of the 
form c(x)m(x) where m(x) is a monomial and c(x) is nonvanishing. Namely, we will use 
Theorem 2.1 to reduce things to proving the following: 

Theorem 3.1: Suppose c(x) is a positive real-analytic function defined on a neighborhood 
of [0, l] n , and m(x) = fllLi X T* ls a nonconstant monomial. Let Z denote {x G [0, l] n : 
m(x) = 0}. Define K^it) by 

K 4> {t)= [ (j)(x)dx (3.1) 

J {xe(0,l) n :c(x)m(x)<t} 

Then if <p(x) is a smooth function on [0, l] n supported in a sufficiently small neighborhood 
of Z, K^it) has an asymptotic expansion of the following form, where the a range over an 
arithmetic progression of positive numbers: 

n-l 

k *v) = EE k *M r in ^y + (3.2) 

There are M > and D a > depending on c(x) and m(x) such that each ki jCe is a 
distribution with respect to 0, supported on Z ', satisfying 

\k ijCe ((f))\ < D a sup sup\dP<f>(x)\ (3.3) 

\/3\<Ma xEZ 

Also, there are e > and C a > depending on c(x) and m(x) such that for any any a and 
any i satisfying < / < a one has 

\^-E a (t)\<C a sup sup \d^(f>(x)\t a+e - 1 (3.4) 

|/3|<M(a+0 x6(0,l)™ 
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Proof that Theorem 3.1 implies Theorem 1.1: Suppose we are in the setting of 
Theorem 1.1. We apply the resolution of singularities algorithm simultaneously to each 
fi(x), each gi(x), and each difference fi(x) — f m (x). Let U be as in the version of the main 
theorem of [Gr] from Theorem 2.1, and let A = {x G R n : gi(x) > 0, ...,gk(x) > 0}. Let 
rj G C C (U) such that r\ = 1 on a neighborhood V of the origin, and let t] = J2iVi be the 
decomposition given by the version of the main theorem given by Theorem 2.1. 

For any (/) G C C (V), observe that = J^. 07^. Thus we have 



J t = I (j)(x) dx = j> I 

J{x£A:0<f 1 (x)<t,...,0<f l (x)<t} i J{ 



4>(x)rji(x) dx 

'{xeA:0</i(i)<t, ... ,0</;(x)<t} 

Let be the composition of coordinate changes corresponding to r\i, let Bi = ^~ A, and 
let Ri be the rectangle (0, a\), (0, a l n ) given by Theorem 2.1. Let Fj = /j o ^j. Then the 
ith term of the above expression becomes 

/ 0o* i (x)(^o* i (a;))J i (a;)da; (3.5) 

^{x6B l ni? i :0<Fi(x)<t, ... ,0<F((x)<t} 

Here Ji(x) denotes the Jacobian and by Theorem 2.1 each r\i o ^(x) is smooth. Since 
each is being resolved by ^i(x) : in (3.5) each ^ o ^i(x) is of the form ci(x)mi(x), 
where mi(x) is a monomial and is nonvanishing. In particular each g\ o ^i(x) is 

either everywhere positive or everywhere negative on Ri. Thus we have that either some 
Bid Ri = 0, whereupon (3.5) is zero, or that each Bid Ri = Ri, whereupon (3.5) is equal 
to 

/ 0o* i (x)(r 7i o* i (x))J i (x)da ; (3.6) 

J {xeR i :0<F 1 (x)<t, ... ,0<F l (x)<t} 

Clearly we only have to consider terms of the form (3.6). Next, notice that since we 
resolved all differences fi(x) — f m (x), each Fi(x) —F m (x) is also of the form ci m (x)mi m (x) , 
where mi m (x) is a monomial and ci m (x) is nonvanishing. As a result, each Fi(x) — F m (x) 
is either everywhere positive or everywhere negative on Ri. Hence there is some Fi(x) 
which is strictly larger than every other F m {x) everywhere on Ri, and (3.6) is 



(p o q^x)^ o q t (x))Ji(x) dx 

{xeRi-.0<Fi(x)<t} 

Since the singularities of gi(x) are also being resolved, we can write Fi(x) = c(x)m(x) and 
the above becomes 

/ (po^ i (x)(r] l o^ l (x))J l (x)dx (3.7) 

J {xeRi:0<c(x)m(x)<t} 

To get (3.7) into the form (3.1), we do a scaling in (3.7) to convert Ri into (0, l) n . Denoting 
this scaling coordinate change by (x±, x n ) — > sx = (sixi, s n x n ), (3.7) becomes 

n „ 

J\si / (p o^i(sx)(r]i o^i(sx))Ji(sx) dx 

l=1 J {xe(0,l) n -0<c(sx)m(sx))<t} 
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This integral is of the form (3.1), with fa replaced by 0o ^ i (sx)(rj i o ^i(sx)) Ji(sx). Hence 
applying Theorem 3.1 will give an asymptotic expansion for (3.7) satisfying (3.2) — (3.4) for 
(f)oty i (sx)(r}ioty i (sx)) Ji(sx) in place of 4>i(x). By the chain and product rules, (3.2) — (3.4) 
then also hold for <p itself. Adding over all i gives Theorem 1.1 and we are done. 

Proof of Theorem 3.1: We induct on the number of positive nii in m(x) = Yli 00 ^- 
We start with the case where exactly one mi is positive, and without generality we may 
assume % = 1, so that c(x)m(x) = c{x)x\ for some k. Since c(x) is a positive real analytic 
function, we can write c(x) = ci(x) k for some positive real analytic ci, and we have 
c(x)m(x) = (ci(x)xi) k . We now do a coordinate change as follows. For i > 1, we let 
Hi = Xi. For i = 1, we let y\ = ci(x)xi. This is a smooth coordinate change on x\ < 5 
for an appropriate 5. Denote the coordinate change map by x = ^(y)- Assuming fax) is 
supported on x\ < 5, we have 

K+(t)= [ 1 faV(y))J(y)dy (3.8) 

J {xE(0,l) n :yi<tk} 

Here J(y) denotes the Jacobian of this coordinate change. Thus K^{t) is the indefinite 
integral of <E>(y) = fa^{y))J{y) in the y\ variable from to t\' k . Hence if we substitute 
faV{y))J{y) = E^o f|(0,2/2,...,yn)ff + 0(y? +1 ) into (3.8), one obtains 

K*(t) = }^( / — T (0,y2,...,y n )dy 2 ...dy n )— — - + 0(t— ) 
fr£ 7(0,1)"-! dy[ (/ + 1)! 

This gives the desired asymptotic expression (3.2) for K^t) in powers of . The ex- 
pressions (3.3) and (3.4) follow from the chain rule applied to <E>(y) = fa^>{y))J{y). This 
concludes the proof of Theorem 3.1 for the case that only one rrii is positive. 

Next, we assume that some / > 1 of the are positive, and we have shown the 
result for / — 1. Without loss of generality, once again we assume that mi = k > 0. Like 
before we let c\{x) be such that c\(x) k = c{x), so that c(x)m(x) = (ci(x)xi) k Y\.i>\ X T* ■ 
Like above, for x\ smaller than some 5 we can do the coordinate change to variables yi, 
where y\ = ci(x)xi and where y^ = Xi for i > 1. Thus if x = \I> (y) , we have 

n 

c(y(y))m(y(y)) = l[yr i (3.9) 

i=l 

Let £(t) be a smooth function on [0, oo) that is supported in [0, 5) : and which is equal to 
1 on a neighborhood of 0. Correspondingly we write fax) = fa(x) + 4>2{x), where 

fa{x) = fax)C(xi), fa{x) = fax)(l - e(xi)) (3.10) 

We correspondingly write K^t) = K\{t) + K^(t), where 

K\(t)= [ fa(x)dx, K%t) = [ fa(x)dx (3.11) 

J {x€(0,l) n :c(x)m(x)<t} J{xe(0,l)":c(i)m(i)<t} 
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Note that K^(t) reduces to the case when / — 1 of the rrij are positive, for the integrand is 
supported on x\ > 5' and thus the x\ factor can be incorporated into the c(x): One does 
a linear variable change in the x\ variable to turn [5', 1] x [0, l] n_1 into [0, l] n and then 
applies the / — 1 case. Hence it suffices to restrict our attention to K^(t). We change to 
the y variables, obtaining 

K\{t)= I M*(y))J(y)dy (3.12) 

J{xe(o,i)-.U" =1 yp<t} 

J(y) is the Jacobian of the coordinate change. In (3.12), the c{x) factor has been removed 
from the domain of integration, a fact that will make our arguments simpler. Write 
= 0i(*(y))J(y), and (3.12) becomes 

Kl(t)= [ <$>{y)dy (3.12') 

Lemma 3.2: To prove K^(t) has an asymptotic expansion satisfying the conditions of 
Theorem (3.2) — (3.4), thereby proving Theorem 3.1, it suffices to find an asymptotic expan- 
sion for K^(t) satisfying the analogues of (3.2) — (3.4) with replaced by $. Furthermore, 
it suffices to consider only the case where rrii > for all i. 

Proof: Since $(y) = (j>i( i S(y))J(y) = £(\I/i(y))0(\I/(y)) J(y), by the chain and product 
rules, |9 a $(x)| is bounded by C a X)|/3|<|a| \d^4>( x )\- Hence if the versions of (3.2) — (3.4) 
hold for K^(t) with $ in place of 0, they will also hold for K^(t) with 4> itself. This gives 
the first statement of Lemma 3.2. As for the second statement, suppose we have proved 
the asymptotics in the case where each rrii > 0. Then in the general situation, we can fix 
those yi variables for which rrii = 0. Then the asymptotics hold for the integral in the 
remaining variables. By integrating the asymptotics with respect to the y, variables for 
which rrii = 0, one sees that the asymptotics hold for the original integral. This completes 
the proof of Lemma 3.2. 

Lemma 3.3: Suppose g(y) is a C°° function on [0, l] n . Let ^ a g a y a denote the Taylor 
expansion for g(y) about the origin. Then for each N we can write 



9(v)= Y 9 a y a + Y q t ^J M^+i'-^n) - Y J?/^ (^yh+2,-,y n ) y ' 



(3.13) 

Here the sum in (3 ranges over all Pi) with < / < n and < < N for each i, lp 

denotes the number of entries f3 has, and hp(yi /3+ i, y n ) denotes <9^y(0, .., 0, y^+i, .., y n )- 

Proof: Taylor expanding in y\ we have 

N-l ~i i N-1 „j j 

9(y) = Y w {0 ^ - Vn)V i + {9{y) ~ £ w (0 ' y2 ' - Vn)V i ] (3 - 14) 

i=0 yi ' i=0 yi 
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This gives the result for n = 1. When n > 1, one substitutes the n — 1 dimensional case 
for |^f-(0,y2, into each term of the left series of (3.14) and the lemma follows. 

Lemma 3.4: Let p(y) denote one of the terms of the second sum of (3.13); that is, let 
p(y) be of the form 

Pi— pi- i=0 c^+i *• 

Let j denote the index called lp + 1 in this term. Then if 7 is a multiindex such that 
7j = 0, then 

|<9M2/)I<C| 7 | ,jv yf\\d\ \c nN +\-*\ 



Proof: If one takes the 7 derivative of p(y), one obtains some terms of the form 

N-l ^ i 

cy m [ q ( yj ,...,y n )-J2 • y n )f] (3.15) 

Here y m is a monomial, and q is some partial derivative of the appropriate hp of order 
at most |7|. Note that the bracketed expression is equal to q(yj, ...,y n ) minus the first 
N terms of its Taylor expansion in the yj direction. As a result, by Taylor's theorem, 

y N a N 

the bracketed expression is equal to ~w~Q^w{Yj, y n ) for some Yj between and yj. 

Hence (3.15) is bounded by c'y^\\q\ \c*. Since q is a derivative of order at most I7I of some 
hp, which is itself a derivative of order at most (n — 1)N of g, we conclude that this term 
is bounded by d'yf\\g\\c^ N +\i\ ■ We add over all terms of d 1 p(y) and we are done. 

We are now in a position to prove Theorem 1.1. 

Proof of Theorem 1.1: By lemma 3.2, it suffices to find an asymptotic expansion for 
K^(t) satisfying the analogue of (3.2) — (3.4) with <p replaced by <E>, and in m(x) = f| i x™ 1 
we may assume m t > for all i. Also, we can assume we know the result inductively for 
dimensions < n. We fix some a > 0. With the goal of finding the asymptotic expansion 
(3.2) up to the power t a , we apply Lemma 3.3 to $>(y), setting equal [(a+e) max^ m^J +1. 
We then insert the result termwise into (3.12'). We obtain two types of terms. The first 
are terms of the form 

f y a dy (3.16) 

Aye(o,D":nr=i^r<*} 

Here $ Q denotes the coefficient of y a in $'s Taylor expansion about the origin. One can 
evaluate (3.16) directly using calculus, by induction on the dimension for example, and 
prove that (3.16) is of the form 

n-l 

$a E E c ^ ln (*r (3.i7) 

p<Ma i=0 
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This is the form that we need for our asymptotics. The second type of term we obtain 
comes from the second sum of (3.13), a term denoted by p(y) in Lemma 3.4. We write 
p(y) = Py ■(&)■> where yj is as in the previous lemma and where y denotes the remaining 
n — 1 variables. We integrate the term first with respect to the y variables and then 
with respect to the yj variable. The term becomes the following, where fh{y) denotes the 



monomial 



/ P yj {y)dydyj (3.18) 

o ^{ye(0,l)™^ 1 :m(j/)<^ r7 } 

We next use the inductive hypothesis on the inner integral, and (3.18) becomes 

,1 n-1 

We expand ln(-Jrj) fc = (ln(t) — rrij \n(yj)) k and (3.19) becomes 

/ [£ E «*to,«(Pv i ){^) a Ht) kl Mvj)*+Ea,„A-^)]dVj (3-20) 

J ° Oi<a fci,fc 2 =0 y j y j 

Here each Kk lt k 2 ,a is a constant multiple of some kk, a - I n particular, each Kk lt k 2 ,a satisfies 
(3.3). We rewrite (3.20) as 

n-1 



^ g ( r 1 K klM AP^)Hy^ dy _ r Ht)kl dyj + r K ^. (321) 

^ 7 , n Jo y^i Jo y^ 

a<a k\,k 2 —0 J J 



We will see that (3.21) gives the desired asymptotic expansion. For by (3.3) and then 
Lemma 3.4 we have 

|«fei,fc 2 ,a(Pj/ J -)l < ^a|bj/jllc M «(Z) 

< D' a yf\\§\\ C M a+n N( Z } (3.22) 
Hence, the absolute value of the t a ln(t) fcl coefficient of (3.21) satisfies 

|> / mjo dyj\ < D a \\$\\ CMa + n N( Z) / ^ ^ — dyj (3.23) 

fe2 

Since N > max/t m^a + e) > rrija > rrija, the right-hand integral above is bounded, so 
(3.23) is at most C' a \\Q\\cM a +n.N ^ z y Since a < a and N < M'a for some M', we have 
Ma + nN < (M + M'n)a, and we conclude that (3.23) is at most 

C'a\\®\\c( M + M ' n ) a (Z) 
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— 1 1^1 lc( M + M '™) a (z) 

Here C' a ' — sup a < a C' a . This gives us the desired estimate (3.3) for the t a ln(£) 1 coefficient. 

We also have to analyze the error term E atVj {^rj) dyj. If < / < a, the Zth 
derivative of this with respect to t is equal to 



d m Eg,y t 



y^-J^(-™j)dyi (3-24) 



Substituting (3.4) into this, this is bounded in absolute value by 



C a fy] ^||^ j || CMa((0)1) „. 1) (4-) a+e - Z (3.25) 
Jo Vj 



By Lemma 3.4 this in turn is bounded by 

^-y f —lm,j-\-N 1 1 -p- I I / t , a _i_ e _; j 

/ Vj \\®\\c Ma + nN {(0,l) n )\^Tj) dVj 

Jo Vj 

= C a m\ CMa+nN{ ^ ir) {f y"- am i- em i d Vj )t a ^- 1 (3.26) 

Jo 

Since A" > arrij + errij, (3.26) is at most 

1 1 ^ I lc Ma +™ JV ((o,i)™)^ a ~'~ e ' (3.27) 
And because A" < M'a for some M', (3.27) is bounded by 

Q(M+M'n)a^(Q 1)™)^ (3.28) 

This gives the desired estimate (3.3) with $ and we are done. To be clear, what we 
showed is that each t a ln(t) J, s coefficient is bounded in terms of C M ° a norms of $ for 
Mq = M + M'n, with corresponding estimates for the error terms. Note that for (3.2)- 
(3.4) to hold we need C M ° a norms, which makes a difference if a is a lot smaller than a. 
However, if for a given a we consider the estimate obtained from the expansion to degree 
a = 2a for example, we get the desired estimates. 
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